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文章 编号 :1005-3085(2011)03-0300-07 


不 定 二 次 规划 全 局 求解 的 一 个 新 算法 * 


IEL, 刘 三 阳 !， 张 建 科 13 
(1- 西安 电子 科技 大 学 理学 院 ， 西 安 710071; 
2- 河南 师范 大 学 数学 与 信息 科学 学 院 ， 新 乡 453007; 3- 西安 邮电 学 院 应 用 数理 系 ， 西 安 710121) 
摘 要 : 针对 工程 设计 、 设 施 布 局 等 领域 出 现 的 不 定 二 次 规划 问题 的 求解 ， 本 文 给 出 了 一 个 新 的 全 局 优化 
算法 . 首先 根据 二 次 函数 的 特点 ， 利 用 线性 松弛 化 技巧 ， 建 立 不 定 二 次 规划 问题 的 松弛 线性 规划 
问题 ， 然 后 通过 一 系列 松弛 线性 规划 问题 的 解 逐步 逼近 原 问 题 的 最 优 解 . 理论 上 证 明了 算法 的 收 
敛 性 ， 数 值 算 例 表明 算法 是 有 效 可 行 的 . 
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本 文 考虑 如 下 不 定 二 次 规划 问题 
(P) min jf(x)-—zTQz-c'zx 
s.t. Ax< b, 
rcX?—-(zeR"|0 <zr<xr<7), 
其 中 Qe R?” 是 nn 阶 实 对 称 和 矩阵，A € R"*"fRIbe R”, cze R. 

二 次 规划 在 工程 设计 、 经 济 学 、 统 计 学 及 设施 布局 等 诸多 领域 有 着 广泛 应 用 Lt,; 站 9， 而 且 许 
多 非 线 性 规划 问题 可 以 转化 为 此 类 问题 的 求解 ， 如 整数 线性 规划 、 整 数 二 次 规划 等 ， 因 此 研 
究 这 类 问题 既 具 有 理论 意义 ， 又 具有 实用 价值 ， 对 于 凸 二 次 规划 的 求解 ， 已 有 很 多 成 熟 的 算 
法 4 对 于 非 凸 二 次 规划 ， 由 于 其 为 N-P 难 问题 ， 所 以 求解 起 来 较为 不 易 . 近 些 年 来 ， 人 们 
对 非 凸 二 次 规划 问题 的 求解 也 做 了 很 多 工作 56 引 ， 本文 利用 二 次 函数 自身 的 结构 特点 ， 构 造 了 
一 个 新 的 分 支 定 界 算法 . 理论 上 证 明了 算法 的 收敛 性 ， 数 值 算 例 表明 算法 是 切实 可 行 的 . 


2 ”线性 松弛 规划 
令 Qi 表示 的 第 i 行 ， 首先 求解 2n 个 线性 规划 问题 


lL min Qiz u;i max Qiz 
s.t. AÁgr<b, s.t. Az <b, 
O<S<TITSTXT, 0zrtizrztzm. 
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并 引进 下 面 的 指标 集 和 符号 
= {i| >0, i=1,2, n}, ed 
Ts = {ilie {L2 VOL JT), z 1n) z-1m(u), iem, 


z; -lni((-w)  zi-ln(-l) i€T» z- S Z; = exp(ui), i€ T3. 


由 于 zi > 0， 所 以 问题 (P) 可 等 价 地 转化 为 如 下 问题 


(EP) min wo(z,z) = 》， z,exp(z;) 十 >》 zi[ zil 一 exp(2i)] 十 5 z,ln(zi) + Tr 
icTi icT» iET3 


s.t. pi(z,z) = Qiz — exp(zi) <0, iem, 
pil£, z) = Qiz + exp(zi) €0, i€ T5, 
pi(z,z)= Qiz — ln(z) <0, ieT,, 

Az <b, 
czeX"-[reR'"|ocgsszzr) 
由 问题 (BP) 和 (P) 的 结构 ， 易 证 如 下 关于 问题 (P) 与 (EP) 的 等 价 性 结论 . 


定理 1 若 (x*,z*) 是 问题 (EP) 的 最 优 解 ， 则 z* 是 问题 (P) 的 最 优 解 ， 反之 ， 若 z* 是 问 
题 (P) 的 最 优 解 ， 则 (x*, z*) 是 问题 (EP) 的 最 优 解 ， 其 中 z* = (2522, sz,» B 


—-In(Qiz", ie Ti, zi -—ln(-Qiz') i€ T5, 2 -—exp(Qiz") iE Ts. 


根据 定理 1， 求 解 问 题 (P) 可 转化 为 其 等 价 问题 (EP) 的 求解 ， 且 问题 (P) 5 (EP) 的 最 优 值 
相等 ， 为 了 求解 等 价 问 题 (EP)， 本 文 给 出 一 个 分 支 定 界 算法 ， 为 此 需要 构造 (EP) 的 松弛 线性 
规划 问题 ， 其 最 优 值 为 问题 (EP) 的 最 优 值 的 下 界 . 

下 面 给 出 (EP) 的 松弛 线性 规划 问题 的 构造 过 程 . 不 失 一 般 性 ， 令 X = [T] CXI 表示 初 
始 和 矩形 或 由 算法 产生 的 子 和 矩形 ， 因 z; > Zi > 0， 所 以 可 构造 问题 (EP) 在 矩形 X 上 的 松弛 线性 
规划 问题 如 下 


(RLP) min wb)(z,z) ENS Zi exp(z;) 十 Xal 2i[ — exp(zi)] + 5 ziln(2;) 十 crz 


i€Ti iET2 icTs 
st. t(x, z) = Qiz — exp(zi) <0, iE, 
p(x, z) = Qiz +exp(z:) <0, i€ T», 
pi(z,z) = Qiz -ln(z)X0, ieTs, 
Az <b, 
reX-(zeR'"|0szzzczz) 
Ti v[(P)] 表示 问题 (P) 的 最 优 值 ， 显 然 有 v[(EP)] > v[(RLP)]. 


定理 2 随 着 lz z zl 一 0， 有 polz, z) 2d quz, z) > 0, pilz, z) n eir, z) 250,i€Ty 
或 ie Ty, Hie Ts. 
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证 明 fi |z — || 一 0 时 ， RE z EX, EAE Z 一 2z; 一 0， 所 以 
polz, 2) — polz, z) = J rilexp(2i) — exp(z:)) + Y, zi(exp(zi) — exp(z:)) 


i€Ti icTs 


十 p» z;(1n(z;) — In(z;)) 0. 
ieT, 
类 似 可 证 ，pi(z,z) — vl(z,2) 90, i€ Ty, 或 ie To, Ri € Ty RAE. 
定理 2 说 明 ， 随 着 分 支 过 程 的 进行 ，wi(z, z) 可 以 任意 逼近 pi(z, z). 


3 ”全 局 优化 算法 


本 节 给 出 求解 问题 (EP) 的 全 局 优化 算法 . 在 算法 中 ， 通 过 求解 一 系列 线性 规划 问 
题 (RLP)， 逐 步 改进 (EP) 的 最 优 目标 函数 值 的 上 下 界 ， 最 终 确 定 出 原 问 题 的 全 局 最 优 解 . 
假定 在 算法 进行 的 第 有 次 迭代 中 ，Gx 表示 由 活动 节点 ( 即 可 能 存在 全 局 最 优 解 的 子 长 方 
TE) 构成 的 集合 . 对 每 一 个 节点 X € Gk， 求 解 线性 规划 问题 RLP (X) 得 最 优 值 ZB(X) = 
V[RLP(X)]. M (EP) 的 全 局 最 优 值 的 下 界 为 了 Bk = min{LB(X),V X € Gr}. 对 任意 的 X € 
Gy, Œ RLP (X) 的 最 优 解 可 构造 出 (EP) 的 一 个 新 的 可 行 解 (构造 方法 参看 算法 )， 进 而 可 检查 
是 否 需 要 修正 (EP) 的 上 界 . 选 定 一 个 具有 最 小 下 界 的 活动 节点 ， 然 后 将 其 分 成 两 部 分 ， 对 每 
个 新 的 节点 求 其 相应 的 解 ， 重 复 这 一 过 程 直到 满足 收敛 条 件 为 止 . 

在 算法 中 ， 对 于 失 形 节点 和 的 剖 分 ， 我 们 采用 对 分 规则 ; 


(G) $ j= argmax{ 歼 一 Zi 一 12 n}; 

(i) 令 T1 = (zz) 

(ii) X! = {z € R" | z; < ti < Ti, i $ j, zj S zj wh X^ = {x € R” |z; <z < 
Ti, i Æj, Yj < zj €x 
通过 剖 分 可 以 得 到 两 个 子 矩形 X1 和 X2， 且 intX1l 站 intX2 20, X! UX =X. 

算法 描述 

步 0 选取 e > 0. 找 出 问题 (RLP) 在 矩形 XX = X? 上 的 最 优 解 z0 和 最 优 值 LB(X0). $ 

z0—1n(Qiz, ie Tj, z2? =ln(-Qix?), ie To, 2z} -—exp(Qia)) ie T5, 
则 (x°, 29) 为 问题 (EP) 的 一 可 行 解 . 置 
LBo= LB(X9), UBo = po(2°,2°). 


如 果 UBo 一 LBo < e， 则 停止 计算 . x° 和 (z0, z0) 分 别 为 是 问题 (P) 和 问题 (EP) 的 e- 最 优 解 . 
fij, &Go-(X9?), F-0, k=1, Sk: 

Fk kl 

步 kl 4UB, = UB. IH EXS 4r 3c UU ME X71 805 A BIA TÓBIE X", XE c 
R^. 4F-FUÜ(X*-!y 

步 k2 DEAE X XU?, MEH SE (RLP) EX = XKt 上 的 最 优 值 LB(X%t) 及 
AUR o (t = 1,2). 令 


ze In(Qiz^5), icT, a ES In(-Q;z*), i €T», zt E exp(Qiz^^), i € T5, 
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则 (z^, zz 为 问题 (EP) 的 一 可 行 解 . 修正 下 界 
U B, = min (UE, polet, 2*t)), 


并 令 (a, zF) 表示 为 满足 UB = qo (z^, z^) 的 点 ; 
步 k3 WR LB(X**) > UB; (t = 1,2), 那么 令 下 = FU(X*5) 
步 k4 4 
F = F|] {X € Gi | LB(X) > UB; 


步 k5 4 
Gi = {X |X e (Gii (6, x*]), x e F). 


步 k6 $ LB, = min(LB(X) | X € Gy), JF X* e Gx 为 满足 LIBk = LB(X*) HTE 
形 . 如 果 UBk — LBy € e， 则 停止 计算 . z* (z^, z^) 分 别 是 问题 (P) 和 问题 (EP) 的 ec 最 优 
解 . 否则 ， 4k = k+ 1, 转 步 k. 


4 ”算法 收敛 性 及 应 用 


定理 3 假设 问题 (EP) 存在 最 优 解 ， 则 算法 或 者 在 有 限 步 内 得 到 问题 (EP) 的 最 优 解 ， 或 者 
产生 一 无 穷 可 行 解 序列 {(x*, zz)}， 其 聚 点 是 问题 (EP) 的 最 优 解 . 

证 明 ” 当 算法 有 限 步 终止 时 ， 结 论 显然 ， 下 证 算法 迁 代 次 数 无 限时 的 收敛 性 . 

根据 文献 [9] 知 ， 一 个 算法 收敛 到 问题 的 全 局 最 优 解 的 充分 条 件 是 界 运算 要 求 满足 
“consistent ”特性 及 界 选取 要 求 满足 “improving” 特 性 . 

所 谓 界 运算 consistent 是 指 在 每 一 步 ， 任 一 未 被 删除 的 部 分 可 被 进一步 剖 分 ， 且 任 一 无 限 被 
剖 分 的 部 分 满足 

„lim (UB, - LBy) = 0, (1) 


XB UB, 是 第 大 次 迭代 在 某 个 子 和 矩形 上 的 上 界 , LB 是 第 有 次 迭代 时 的 最 好 下 界 , LBk, UB 不 
必 同 时 出 现在 同一 子 矩形 上 ， 下面 说 明 (1) 成 立 . 

因为 算法 中 采用 的 矩形 对 分 是 穷 举 的 ， 所 以 根据 文献 [9] 知 ，(1) 成立， 即 算 法 中 的 界 运算 
是 consistent. 

界 选取 满足 improving 特性 是 指 在 有 限 次 剖 分 后 ， 至 少 有 一 个 下 界 在 其 上 达到 的 矩形 被 选 
出 ， 作 为 进一步 剖 分 的 矩形 ,根据 算法 ， 在 迭代 中 ， 作 为 进一步 划分 的 矩形 恰恰 是 下 界 在 其 上 
达到 的 矩形 ， 因 此 界 选 取 是 improving 的 . 

综 上 可 知 ， 本 文 给 出 的 算法 满足 界 运 算是 consistent 的 ， 且 界 选取 是 improving 的 ， 因 此 根 
据 文献 [9] 中 定理 IV.3. 知 ， 该 算法 是 全 局 收敛 的 . 

为 了 验证 本 文 方法 的 可 行 性 ， 我 们 利用 Matlab7.1 编 程 ， 做 了 一 些 数 值 实验 . We = 
1.0e 一 3， 变量 个 数 记 为 n， 算 法 运行 时 间 以 秒 为 单位 . 

例 1 


min z7Qz--c?z 
s.t. Art€b, 


z € X? = (0.4743 € xı € 5.8649, 0 < z2 < 5.0279,0 € za < 2.5783), 
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其 中 
c = (—0.992372 — 0.046466 0.891766)7, 
b = (2.865062 — 1.491608 0.519588 1.584087 2.198036 — 1.301853 — 0.738290), 


0.992934 —0.640117 0.337286 
Q = | —0.640117 —0.814622 0.960807 |, 
0.337286 0.960807 0.500874 


0.488509 0.063565 0.945686 
—0.578592 —0.324014 —0.501754 
—0.719203 0.099562 0.445225 
A= | —0.346896 0.637939  —0.257623 
—0.202821 0.647361 0.920135 
—0.983091 —0.886420 —0.802444 
—0.305441 —0.180123 —0.515399 
最 优 解 为 (5.2107 5.0279 0.0000), RIA -32.5794, R1, 运行 时 间 0.401768 fb. 
5/2 
min z'/Qz- cz 
s.t. Arctb, 


zeX?-(0€z21€2,0€z5«2), 


其 中 

-4 2 
0 1 

-1 2 
ce-(247, b=(1 2 43 1, Q= , A=| 1 1 

2 —4 
1 0 
1 -4 


最 优 解 为 (0.75 2), S&ULIB 7 -1.0625, W3 次 ， 运 行 时 间 0.310939 $. 
下 面 我 们 对 两 个 大 规模 的 标准 非 正定 二 次 规划 问题 进行 测试 . 
例 3 


n 
: 1 
min FX (3? T zm) 
i=1 
s.t. z EX? = {0< ri <1, i=1,2, n}. 


例 3 是 一 个 凹 可 分 二 次 规划 问题 ， 表 1 给 出 了 例 3 的 测试 结果 . 


305 




















第 3 期 汪 春 峰 ， 等 : 不 定 二 次 规划 全 局 求解 的 一 个 新 算法 
表 1: 例 3 的 测试 结果 

n RAA 运行 时 间 最 优 值 
30 1 0.852020 -45 

50 1 1.628337 -75 
80 1 3.098150 -120 
100 1 4.270962 -150 
150 1 8.128278 -225 
200 1 13.375105 -300 

例 4 
n 
min — > T; 
i=1 
j 
s.t Mus j212 n, 
i=1 
z;> 0, i=1,2, 
例 4 的 测试 结果 如 表 2 所 示 . 
A2: 例 4 的 测试 结果 

n 和 迭代 次 数 运行 时 间 最 优 值 
30 1 1.320828 -900 
50 1 3.592944 -2500 
80 1 10.694849 -6400 
100 1 20.450738 -10000 
150 1 73.106514 -22500 
200 1 193.654163 -40000 








由 以 上 数值 算 例 可 见 ， 我 们 的 算法 是 有 效 可 行 的 . 
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A New Global Optimization Algorithm for Indefinite 
Quadratic Programs 


WANG Chun-feng!?, LIU San-yang!, ZHANG Jian-ke!^? 


(1- School of Mathematical Sciences, Xidian University, Xi'an 710071; 
2- Department of Mathematics, Henan Norma] University, Xinxiang 453007; 3- Department of 
Mathematics and Physics, Xi'an Institute of Posts and Telecommunications, Xi'an 710121) 


Abstract: A new global optimization algorithm is presented to solve indefinite quadratic programming 
problems which have been extensively used in engineering design and facilities layout etc. First, the 
linear relaxation programming problem of the indefinite quadratic programming problem is derived by 
utilizing the characteristics of the quadratic function and a linear relaxation technique; then, by means 
of the sequential solutions of a series of linear programming problems, the global optimal solution is 
obtained. The theoretical analysis and numerical experiment show that the presented algorithm is 
convergent and efficient. 

Keywords: indefinite quadratic programming; global optimization; linearization relaxation; branch 
and bound 
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